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ABSTRACT 

HETE-II and BeppoSAX have produced a sample of GRBs and XRFs with 
known redshifts and Epk- This sample provides four important empirical con- 
straints on the nature of the source jets: Log Eiso is approximately uniformly 
distributed over several orders of magnitude; the inferred prompt energy Log E^ 
is narrowly distributed; the Amati relation holds between Ei^o and Ep^] and the 
Ghirlanda relation holds between E^ and Ep^. 

We explore the implications of these constraints for GRB jet structure during 
the prompt emission phase. We infer the underlying angular profiles from the 
first two of the above constraints assuming all jets have the same profile and total 
energy, and show that such "universal jet" models cannot satisfy both constraints. 

We introduce a general and efficient method for calculating relativistic emis- 
sion distributions and Epk distributions from jets with arbitrary (smooth) angular 
jet profiles. We also exhibit explicit analytical formulas for emission from top-hat 
jets (which are not smooth). We use these methods to exhibit Ep^ and Eiso as a 
function of viewing angle, for several interesting families of GRB jet profiles. We 
use the same methods to calculate expected frequency distributions of Eiso and 
E^ for the same families of models. 

We then proceed to explore the behavior of universal jet models under a range 
of profile shapes and parameters, to map the extent to which these models can 
conform to the above four empirical constraints. 

Subject headings: Gamma Rays: Bursts — ISM: Jets and Outfiows — Shock 
Waves 



1. Introduction 



Since the recognition that gamma-ray burst (GRB) sources are highly relativistic coUi- 
mated jets, the study of the structure of the jets and of the distribution of jet parameters 
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among burst sources has been an area of active research. An early observational success 
of the jet picture of GRBs was the discovery by Frail et al. (2001) of a striking correlation 
between GRB energy fiuence and the break times of GRB afterglow light curves. Using the 
simple picture of a conical, uniform jet, Frail et al. (2001), and later Bloom et al. (2003) 
inferred jet opening angles from the break times using results on the dynamics of after- 
glows of Sari et al. (1999) for a sample of GRBs with known redshifts. After correcting the 
isotropic-equivalent gamma-ray energy of each GRB by the inferred solid angle subtended by 
its conical jet, Frail et al. (2001) found that the extremely broad (3-4 decades) distribution 
of isotropic-equivalent energy was mapped to a remarkably narrow (~ 1 decade) distribution 
of inferred total energy. 

The debate over the meaning of this result has given rise to two main alternative inter- 
pretations. The first interpretation is an elaboration of that of Frail et al. (2001), in which 
GRB sources all have a common amount of energy available for gamma-ray emission (they 
are "standard bombs"), and differ from one another by the opening angle of their conical 
jet (Lamb, Donaghy, & Graziani 2005). Implicit in this view is that little or no emission is 
received by an observer whose line-of-sight is not included in the jet cone. 

The alternative view is that GRB jets might have structure that is richer than a simple 
uniformly-emitting cone. The attractive feature of this picture is that if the jet structure is 
sufficiently rich, it might be possible to ascribe the same jet structure to all GRB jets, and 
to explain the distribution of GRB energies using only the variation of observer lines-of-sight 
with respect to the jet axis (Rossi et al. 2002; Zhang & Meszaros 2002). 

These two viewpoints are sometimes presented juxtaposed as the "uniform jet model" 
and the "structured jet model". This terminology is unfortunate, since it is clear that the 
real distinction under debate is not really whether GRB jets have "structure", but rather 
whether the observed energetics of GRBs can be explained in terms of a single parameter (the 
angle between the observer's line-of-sight and the jet axis), or whether additional parameters 
relating to jet structure are required. Thus the juxtaposition would be better thought of as 
the "universal jet model" versus the "variable-geometry jet model," where a sort of middle 
ground is occupied by the "quasi-universal model" ( Lloyd- Ronning et al. 2004). 

One element that has not featured prominently in discussions of jet structure is the 
role of relativistic kinematics in prompt GRB emission. Relativistic kinematics is a central 
feature of the "off-axis beaming" models of Yamazaki and collaborators (Yamazaki et al. 
2002, 2003, 2004; Toma et al. 2005), which those authors used to attempt to unify the 
phenomenon of GRBs with that of "X-Ray Flashes" or XRFs (Heise et al. 2001; Kippen 
et al. 2001). However, their models invoke a uniform conical "top-hat" jet structure. By 
and large, calculations of prompt emission expected from "structured" jets have essentially 
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bypassed the role of relativistic beaming and Doppler shift, assuming a "what you see is 
what you get" relation of the jet emission profile to the jet radiation pattern that tacitly 
implies a Lorentz factor satisfying 7"^ (Zhang & Meszaros 2002; Lloyd-Ronning et al. 
2004; Zhang et al. 2004). 

This is a gap that is worth plugging. It is clear that the observed radiation pattern 
from a "structured" GRB jet can be substantially modified by relativistic kinematics, both 
spectrally and in terms of total radiated energy. Depending on the Lorentz factor and viewing 
angle, relativistic kinematics can turn a steep jet profile into a shallow radiation pattern, or 
turn a hard source spectrum into a soft observed spectrum (or, as we shall see, vice versa). 
It is certainly not safe to assume that it is even possible to produce a given radiation pattern 
without a stringent — and possibly unphysical — constraints on the Lorentz factor. 

Certain constraints on universal jet shapes and distributions may be inferred from ob- 
servations. The empirical constraints on jet models that we consider here are illustrated in 
Fig. 1. They are: 

• The distribution of GRB isotropic-equivalent energy Ei^o appears to be quite broad, 
spanning at least four orders of magnitude, and appears approximately uniform in 
logE'iso, that is, 

dN/d\ogEi,^^A, (1) 
where A is a constant (Amati et al. 2002; Lamb et al. 2004). 

• The distribution of GRB total gamma-ray energy inferred using the Frail et al. (2001) 
procedure — E^^ — appears to be quite narrow, with most values clustered within a 
decade around a geometric mean whose (model-dependent) value is somewhere in the 
range 5 x 10™-3 x 10^^ erg (Frail et al. 2001). 

• Empirically, GRB isotropic energies Eig^ and peak energies -Epeak of GRB z/F^ spectra 
appear to be satisfy the "Amati Relation" , a very tight correlation according to which 
^peak eIJ^ (Amati et al. 2002; Lamb et al. 2004). 

• Empirically, GRB inferred total energies E!p and peak energies -Epeak of GRB uFi, spec- 
tra appear to be satisfy the "Ghirlanda Relation" , a very tight correlation according 
to which £;peak ~ -^7^° (Ghirlanda et al. 2004). 

In this paper we consider the extent to which universal jet models with various under- 
lying emission profiles can satisfy these observational constraints on GRBs with afterglows 
and redshifts. The plan of the paper is as follows: 
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In §2 we consider the effect of the choice of functional form of the emission profile on the 
frequency distributions of E^so and We consider power-law profiles, and introduce the 
"Fisher" (exponential) profile, which is designed to satisfy Eq. (1) exactly. We also inquire 
into the distributional properties of top-hat jet profiles in universal jet models. 

The treatment in §2 uses "bare" emission, that is, uncorrected by relativistic kinemat- 
ics. In §3 we introduce a general and efficient method of computing fully relativistic prompt 
emission from jets with arbitrary (smooth) emission profiles. We introduce the notion of 
"photon-number- weighted average Doppler shift", a quantity that serves as a proxy for ob- 
served Epcak- We also exhibit the method for calculating how the frequency distributions of 
§2 are modified by relativistic kinematics. Finally, we give compact analytic expressions for 
relativistically-correct emission and average Doppler factor in the special case of the top-hat 
jet profile. 

In §4, we use these ideas to explore the characteristic angular emission profiles, Eiso-Epeak 
("Amati"), and E^^-Epeak ( "Ghirlanda" ) relations that would be observed from Universal 
GRB models with Fisher, power-law, and top-hat emission profiles. We also show exphcitly 
how the frequency distributions of §2 are modified by relativistic kinematics. 

Our conclusions concerning the ability of universal jets based on the various emission 
profiles to satisfy the observational constraints are discussed in §5 

2. Universal Jet Profiles and Their Frequency Distributions in Energy 

For the sake of clarity of the discussion that follows, it is useful to draw a terminological 
distinction between jet profiles and jet models. 

A jet profile is a distribution of emission across the surface of the emitting shock. Popular 
choices include the "top hat" jet, and 9^^ power-laws of Rossi et al. (2002) and Zhang & 
Meszaros (2002) (also sometimes referred to as the "structured jet"), where 9 denotes angular 
location on the jet. 

A jet model of GRB emission includes a choice of profile, but also a choice of distribution 
of profile parameters. For example, the model of Zhang & Meszaros (2002) selects a power- 
law with a fixed index, and a fixed total energy. The "Variable Opening Angle" model of 
Lamb, Donaghy, & Graziani (2005), on the other hand, selects a top-hat profile and assumes 
a fixed total energy, but a distribution of jet opening angles. 

Models such as that of Rossi et al. (2002) and Zhang & Meszaros (2002) that adopt 
a fixed profile shape and total prompt energy — and thus rely only on the variation of 
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observer viewing angle to produce all the distributions of GRB energetics — are properly 
called "universal" models. We will refer to models that allow for variations in profile shape 
as "variable-geometry" models. Note that neither class of models necessarily requires that 
its underlying profile should belong to a particular family — one can in principle "mix-and- 
match" models and profiles quite freely. 

It is apparent that universal-jet models of GRBs are simpler than variable-geometry 
models, in the sense that they require fewer parameters. Clearly, an important question that 
must be addressed in this connection is: do universal jet models have enough freedom to 
satisfy known observational constraints on GRBs? This question is the main concern of this 
paper. 

As a notational preliminary, let Xiso = lni?iso) and yU^, = (1 — cos^^^), where 6^ is the 
angle between the line-of-sight and the center of the jet. We will also use = InE^ 

and x!^'' = lnE^\ where "(T)" refers to "True" (i.e. the actual total gamma-ray energy 
emitted), and "(I)" refers to "Inferred" — the gamma-ray energy inferred using the Prail et 
al. (2001) procedure, in which the jet energy is corrected from its "isotropic" value using the 
solid angle subtended by the jet opening angle inferred from the break time of the jet. In 
general, E? ^ e!^^\ 

In this section, we assume symmetric "2-lobed" jets, so that < /ly < 1. 



2.1. The Fisher Jet Profile 

Suppose we describe a Universal Jet by a function f{iJ,v), that is 

Xiso^fM- (2) 

How do we choose /(//^) so as to satisfy the constraint of Eq. (1)? 

By construction. Universal jet models assume all properties of GRB energetics are pro- 
duced by variation of the viewing angle. Naturally, this variation is assumed uniform in /ly 
ior < /ly < 1. In this range of we therefore have 

dN _ dN/di^y 
dXiso dXiso/ dfj,y 

= [/'(/..)]-' 1,,,.,.,^. (3) 

where B is an uninteresting normalization constant. This kind of expression corresponds 
closely to expressions of Guetta et al. (2005) concerning universal jet models, the obvious 
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difference being that here we concern ourselves with modehng the E'iso distribution, and not 
the observed fluence distribution. 

In consequence of the hmited allowed range of /ly, dN/dX^so can only be constant in a 
limited range of Xsso- Within that range, we have = const., so that 

^iso = filJ'v) = Cifiv + Co, (4) 

and 

-— = -e(Xi,o-co)0(co + ci-Xi,o). (5) 

The shape of the jet in hnear energy space is 

E^iso = exp(ci/i^ + Co), (6) 

which peaks at fiy = 1 for Ci > 0. The larger Ci, the sharper the peak. This functional form is 
called the "Fisher Distribution" , and is well known from the theory of statistical distributions 
on spheres (Mardia 1972, p. 228). In the small-angle approximation, when /i„ ^ (1 — 6'„^/2), it 
approaches a symmetric 2-dimensional Gaussian with standard deviation ^o = ci~^/^, a form 
considered by Zhang & Meszaros (2002), Zhang et al. (2004), Lloyd- Ronning et al. (2004), 
and Dai & Zhang (2005). The Fisher distribution is more convenient for our purposes, since 
its geometrically-natural form makes it useful for analytical work for widths and viewing 
angles well outside the small-angle approximation. 

Note that the domain restriction in Eq. (5) says that Xiso varies over a domain whose 
breadth is ci = l/9o^. For a jet with an opening angle of 5.7°, 6^0 = 0.1, and the resulting 
range of i^iso is 100 e-foldings, corresponding to about 43 decades. Naturally, instrument 
thresholds prevent us from probing such a range of £^iso- 

It is convenient to adjust the normalization of this expression so that it reflects the true 
emitted energy, e!^\ Since E'iso = "indE^^ /d'^Qy, we have 

Ef) = 2 f -dii,e''^l''"+'^ 
Jo 2 

= e^°^o' (exp(0 - 1) , (7) 

where the factor of 2 in the first line accounts for the energy in both lobes of the symmetric 
jet. 



We thus have 



9q (exp(^o ) - 1) 



- 7- 



and 




(9) 



2.2. The X!^ ' Distribution of Universal Fisher Jets 



The quantity E^^o is directly related to a physical property of the GRB — the energy 
emitted per unit solid angle along the observer's line-of-sight. The quantity E^^\ on the 
other hand, has a shghtly different status. It is the result of a data-analysis procedure — the 
Frail et al. (2001) procedure — which multiplies Eiso by a solid angle interior to a jet opening 
angle inferred from the light-curve break time. This procedure is particularly adapted to 
top-hat profile jets, so the physical meaning of e'^^ in the context of other profiles can 
be somewhat contrived. Nonetheless, the empirical distribution of e'^^ appears to supply 
important constraints on jet models, as noted above, so it is important to understand jet 
model predictions of these distributions. 

In the case of top-hat jet profiles, the Frail ct al. (2001) procedure infers a jet opening 
angle from the afterglow light-curve break time. In the case of a non-uniform universal jet, 
the break time docs not correspond to a jet edge. Rather, it occurs when the bulk of the 
emission from the jet is accessible to the observer, since after that time, the evolution of the 
light curve is controlled by the evolution of 7 only, which is expected to be a power-law in 
time. 

As pointed out by Kumar & Granot (2003) and Dai & Zhang (2005), the light-curve 
break indicates one of two situations: (1) the viewing angle 6^ is larger than 6q, and the 
beaming angle I/7 has widened to the point that the peak of the jet has become visible, 
so that the emission can no longer increase in brightness at the rate that it did when the 
beaming cone was still climbing up the peak. Or (2) the viewing angle 9^ is smaller than 
^0, and the beaming angle I/7 has widened until it is equal to ^0, so that most of the jet 
emission is now visible, and further increase in 1 produces a diminishing return in flux. 

Following Kumar & Granot (2003) and Dai & Zhang (2005), we therefore assume that 
the inferred break angle 9}yr — max(^^, ^0), and that E^^ is given by 



Em=Ei,oX(l-/X;.,), 



(10) 



where ji^yr = cos 
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We may now inquire as to the distribution dN/ dX^^ of Universal Fisher jets. X^^ is 
given by 

^^'^=^iso + ln(l-//6.). (11) 

The shape of this function is shown in Fig. 2. From the figure, it is apparent that the range 
of E^^ is about the same as the range of E^soi about exp(^o~^)- It also appears that the 
correspondence from x'^'^ to ji^ is single-valued in almost all of this range, except for an 
small annulus of inner radius Oq and outer radius -\/2^o, where it is three- valued. 

We may therefore write 

dN dN/dn^ 



E 



= bxj: 

Roots 



1- Hv 



(12) 



where the sum is over the roots of Eq. (11), of which there is only one outside the aforemen- 
tioned annulus. 

Assuming a typical 9q < 0.1, it is apparent that in most of the range < /ly < 1, 
the first term in the summand of Eq. (12) dominates the second term, so that dN/dX^^ is 
approximately also uniform in most of its range. It does spike (integrably) when 1 — ji^ — Oq' 
(as 9y crosses the inner radius of the annulus), but this behavior is not enough to make the 
distribution a narrow one. Most of the probability, and therefore most of the events, occur 
at a range of angles considerably larger than the core annulus, which subtends a small solid 
angle. For example, the 95% point of the dN/dXj^ distribution occurs at a value of Xj'^ 
corresponding to jj,^ = 0.05. With 6q^^ = 100, this value of Xj^ is lower than the value 
near the jet core by about 90 e-foldings, equivalent to about 39 decades. Examples of the 
behavior of Eq. (12) are shown by the "bare" cuves in the four panels of Fig. 9, for a range 
of values of ^o- 

We see that Eq. (12) represents an extremely broad distribution, far broader than can 
actually be probed by real instruments - obviously, HETE, BeppoSAX, BATSE and SWIFT 
didn't and won't probe 40 decades of Ei^^or Eiso- The key point, however, is that in the 
observationally accessible region, the distribution is broad, not narrow, and is only cut off 
by instrumental effects. 
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2.3. The Xiso Distribution of "Universal Structured" Jets 

The universal "structured" jet model of Rossi et al. (2002) and Zhang & Meszaros (2002) 
does not attempt to satisfy constraint (1). Instead, it attempts to preserve the appearance 
of the Frail et al. (2001) result — that is, to satisfy constraint 

dN/dX^^^ = iVo X 5(XW - XW). (13) 

by imposing the following shape: 

^iso = i^/(l-/x.), (14) 

or, equivalently, 

Xi,o = Q-ln(l -//,), (15) 

where K and Q are constants, and where we have replaced the small-angle form with the 
general, spherically-correct form. 

Since this form does not attempt to constrain dN/dX^so, it is of interest to calculate 
what this distribution might be. This is easily done: 



dN _ dN/df.^ xe(x,o-g) 



dXiso dX\so/ djiy 

= s X (1 - //^) X e(Xiso - Q) 

= sKe-^iso X e(Xi3o-g). (16) 

In other words, the distribution is exponential in Xiso, with a width scale of one e- 
folding (0.43 decades) and a cutoff at some low value of Xiso- Needless to say, the empirical 
distribution (Fig. 1) looks nothing like this. 



2.4. Jet Distributions, Broad and Narrow 

It appears from the foregoing that universal jet models have a serious problem: if you 
try to make a jet that has a broad £^iso distribution, you find that the jet also has a broad e'^^ 
distribution, while forcing a narrow E'^'' distribution results in a narrow E^^o distribution. 

To illustrate the point, suppose we let the functional relation between E'iso and l-//„ be 
a general power-law, so that 

^iso = Q-(^ln(l-/x„). (17) 
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The universal "structured" jet model of Rossi et al. (2002) and Zhang & Meszaros (2002) 
corresponds to 5 = 1. Those authors also considered general power-law jet profiles as well, 
although those profiles were power- laws in 9y, rather than in the solid angle interior to 9^ as 
here. 

It then follows, by the same argument that led to Eq. (16), that the energy distributions 

are: 

^^Cx r^e-^iso/^ X e(Xi,o - Q), (18) 

and 

^ = C X (5 - ^ Q(^(I) _ (^9) 

Note that for the Zhang & Meszaros (2002) universal jet profile, 5 — > 1, and the right- 
hand side of Eq. (19) tends to a 5-function. 

It appears from Eqs. (18) and (19) that the family of models satisfying Eq. (17) is 
capable of producing energy distributions of any breadth, simply by varying the value of 6. 
The larger the value of S, the broader the distribution. This property is in contrast with the 
properties of Fisher jets, which (by design) can only represent broad E'iso distributions. We 
may therefore use the power-law jet family of models to explore the connection between the 
breadth of the Xiso distribution and that of the X^'' distribution. 

Prom Eqs. (18) and (19) we may observe that within this family of models, the widths of 
the two distributions are closely linked: if we try to push out the width of the Xig^ distribution 
by increasing 5, we will also automatically increase the width of the x!^^ distribution to very 
nearly the same width. If, on the other hand, wc attempt to restrict the range of X^^ by 
making S close to 1, we thereby also restrict the range of Xiso- 

It appears, then, that Universal jet models force us to choose between two observational 
constraints on GRBs. Within the context of a universal jet model, it is not possible to satisfy 
both the requirement that dN/dXi^o should be broad and the requirement that dN/dX^^ 
should be narrow. 



3. Computing Relativistic Emission PVom Jet Profiles 

In this section we present formulas and methods of computing various observational 
quantities given certain jet angular emission profiles. Naturally, underlying emission profiles 
are not observed directly, but rather they are modulated by relativistic effects due to their 
relativistic bulk motion. 
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In this work, we are concerned with jet emission profiles that are due to variations in 
emissivity across the jet, rather than due to variations in the Lorentz factor 7. We consider 
that during the prompt emission phase of the GRB, 7 is constant across the jet and has no 
appreciable time variation (Rees & Meszaros 1994). 

We consider prompt energies and observed fluences, rather than prompt peak luminosi- 
ties and observed peak fluxes. Given the assumed invariance of 7, this approach allows us to 
ignore time variations in the emission, and to represent prompt emission properties of GRBs 
simply in terms of weighted integrals of jet angular profiles. 

We also assume that the observer's detection process is approximately bolometric — that 
is, we ignore the spectral effects of rclativistic kinematics, so we do not concern ourselves with 
the kinematic transformations of detector bandpasses, or with the integration of emission 
spectra over those bandpasses. 

Note that in this section, unlike in §2, we assume "single-lobed" jets, for the sake 
of simplicity of our expressions. A symmetric "two-lobcd" jet may be constructed from 
oppositely directed single-lobed jets, ascribing half of the total energy to each lobe. 



It is obviously of interest to calculate the observed fluence as a function of viewing angle, 
from which ii^iso follows directly. 

In what follows, we refer to the inertial frame at rest with respect to an element of the 
jet as the "rest frame" (of that element), and to the inertial frame at rest with respect to 
the central engine as the "source frame." 

Suppose a jet with bulk Lorentz factor 7 radiates a soiucc- frame energy E-^ . Let 
each differential element of the jet d'^n be labeled by the unit vector n, so that n furnishes 
coordinates on the jet surface. For concrctcncss, we will assume n is a source-frame (as 
opposed to rest-frame) vector. Define a source-frame unit vector / that points in the direction 
of the line-of-sight to the observer. 

The element of the jet at n radiates isotropically in its own rest frame an amount of 
rest-frame energy per unit rest-frame solid angle 



3.1. 



Observables 



3.1.1. Eiso and K 




An 



X e{n ■ b) d^fi, 



(20) 
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Here is a unit direction vector pointing in the direction of the hne-of-sight to the observer 

— * 

in the the rest frame of the jet element at ft. The profile function e{n-b), which is azimuthally 
symmetric about the direction of the jet-axis unit vector b, is normalized so that 



^ e{n ■ b)d:^n = 27: J dxe{x) = l. (21) 

The usual relativistic Doppler shift may be applied to yield the source-frame energy per 
unit source-frame solid angle: 

= i^^x 7-3(1-/5^-0-3^2^ 

= i e{n ■ b) X - f3n ■ lyU^n. (22) 

Since this is an energy, rather than a luminosity, it should be converted to an observed 
fluence by means of a multiplication by (1 -|- z)dL{z)'~^, where ^^(z) is the usual luminosity 
distance. 

The isotropic-equivalent energy is then given by 

dE 
dn 



Ei^ = 47r X 



= ) ^ e{n ■ b) X 7-^(1 - I3n ■ l)-^ d^n. (23) 

We will also be interested in the inferred total energy, E^^ , which is given by 

dE 

= 27r(l-cos^ft,) X ^ 

^ dn 

= ^(l-cos^b.) x£;i,o, (24) 

where — max(^„,^o)) ^-nd 6^ is the viewing angle — the angular distance between the 
observer and the axis of the jet. Again, this assumes all the inferred energy is ascribed to a 
single jet lobe. 



3.1.2. Peak Energy 



Another interesting calculation is the effective E^^ak observed along the viewing direc- 
tion, assuming a common rest-frame E'peak for the entire jet. For universal jet models, this 
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quantity can be placed in relation to Eiso and E^^ , in order to examine whether the empirical 
Epeak-Eiso relation of Amati et al. (2002), or the empirical Epeak-E^^ relation of Ghirlanda 
et al. (2004) can be reproduced. 

Technically, £^peak is the peak of the i/Fi, spectrum obtained by superposing the appro- 
priately Doppler-shifted and suitably Doppler-weighted spectra received from various parts 
of the jet. Performing this superposition explicitly presents difficulties, such as choosing a 
detailed shape of the underlying emission spectrum in such a way as to yield a Band (GRB) 
function (Band et al. 1993) spectrum for the superposition. This is the approach chosen by 
Yamazaki et al. (2002, 2003, 2004). 

Our approach is to calculate the "photon- number averaged i?pcak-" We ascribe a com- 
mon rest-frame E^^^^ to the entire jet, and calculate the average Doppler-shifted -E'pcak in 
the viewing direction, where the average is weighted by photon number fiuence. This is 
tantamount to computing the average of the Doppler shift D = 7~^(1 — (3n ■ according 
to the weighting function e(n • b) x 7~^(1 — j3n ■ /)~^ and applying it to the rest-frame Ep^ak. 

In other words, we want to be able to calculate 

^^pea^ = (/^) X (25) 

where E^^^ is the rest-frame peak energy, and 

_ §<Pne{n-h) x 7-^(1 - /jn ■ ^^^^ 
§ d?ne{n ■ h) x 7"^(1 - (5n ■ 

The accuracy of Eq. (25) as an estimator of the peak of the effective spectrum depends 
on the nature of the underlying rest-frame spectrum. In general, there will certainly be 
systematic uncertainties and biases in {D)E^^^^ as an estimator of observed peak energy. 
Nonetheless, it seems reasonable to expect that in its functional dependence on ^^.y, in its 
correlations with i^iso and -^7^'', and in its frequency distribution in universal jet models, the 
behavior of {D)E^^^^ should not be too dissimilar from the behavior of the peak energy of 
the true composite spectrum. Since these are the properties that concern us, we adopt (D) 
as a proxy for £^peak in what follows. 

One useful special case of (D) is e(n ■ 6) = 1 — isotropic outflow. In this case, it is 
straightforward to show that (D) = 7. The reason this is useful is that this is also the limit 
to which {D) should tend when the profile e varies weakly on angular scales comparable to 
7~^ on the part of the jet that is moving parallel to the hne-of-sight. It thus provides a 
valuable sanity check in practical computations. 
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3.1.3. Relativistic Corrections to Frequency Distributions of Eiso and E^K 

It is interesting and important to understand the extent to which our §2 results on Eiso 
and Ej^ distributions of universal jet models are affected by relativistic corrections. 

What is required is the form of dN/ dEiso and dN/ dE^^ , where instead of being given 
directly in terms of the emission profile as we assumed above, Eis^ and Ej'^ are given by 
Eqs. (23) and (24). 

Since in universal jet models the distribution of Eiso is entirely due to a uniform distri- 
bution in = cos^^, we may write 

dN dN/dfiy const. 



dEiso I dEiso /djj,v\ I o?£^iso / dji^ \ 



(27) 



We therefore need to calculate dEiso/ d^y starting from Eq. (23). This is accomphshed by 
differentiating the Doppler factor in the integrand with respect to /i^. By adopting spherical- 

— * — * 

polar coordinates {O^, 0^) for I with the 2;-axis aligned with b, and differentiating the quantity 

— * 

n ■ I with respect to /ly while holding the 0^ and n fixed, it may be shown that 
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d{n ■ I) n ■ {b — fij) 
diXy 1 - A*^^ 

It follows then that if we differentiate Eq. (23) with respect to /i^ we obtain 

d/.. -/(l-/..2)/^^(l_/3n-/>^ ^ ' 

We also require 



dN dN/djjiy const. 



(30) 



dE'^'^ I dE^^'^ /dfxyl I dE? /dfiy]' 
In Eq. (30) the factor dE^/dpiy may be related to dEiso/ d^y using Eq. (24): 

'^ = \l_l^^)^_\Eiso. (31) 

djJLy 2 djJLy 2 

A special limiting case is the "needle" jet limit, where a narrow jet is observed off- axis 
(Yamazaki et al. 2002, 2003, 2004). In this case, t{n ■ b) is so narrow that it acts as a 
5-function. By Eq. (23), the radiation pattern is then proportional to (1 — /?/x^)~^. It is 
straightforward to show that in this limit, we have 

lim ^Eiso-^'^. (32) 

dlog£;iso 
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Furthermore, using Eq. (24), and using P fa 1, we can also show that 

a log ' 

These limits furnish useful sanity checks, and are interesting in their own right in the context 
of the off-axis beaming model (Yamazaki et al. 2002, 2003, 2004). 



3.2. Spherical Convolution 



We present in this section an efficient and general technique for computing the energy 
emission distributions of GRB jets with arbitrary (smooth) emission profiles. 



3.2.1. Method 

The quantities of interest described in the previous section are all expressed in terms 
of convolutions on the 2-sphere of two functions that are azimuthally symmetric but not 
mutually concentric. That is, they are of the form 

^ (fnf^\n-q^'^) x f^^\n ■ q^^^). (34) 

Now suppose we represent each of the two functions /*^^) and /^^^ on the two-sphere 
by their expansion in spherical harmonics. Since they are azimuthally symmetric, only the 
m = spherical harmonics appear, that is 

/«(rz • g «) = £ /«^PK^ • ? ^^^) (35) 

l=Q 

where the Pi{x) are the usual Legendre polynomials, normalized so that J^^Pi{xYdx — 
2/{2l + 1). By the orthogonality of the Pi{x), it follows immediately that 

= j\xPi{x) x/«(x). (36) 

If we let the ^-direction be along q^^"^ and define spherical polar coordinates (^, 0) for n, 
so that fi ■ q^^^ = cos 6, then in Eq. (35) we may replace Pi{n ■ q^^^) by v^47r/(2Z + l)yj°(^, 0) 
in the expression for f^^\ 
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If we introduce spherical polar coordinates {02, ^2) for the direction vector ' (so that 
cos ^2 = Q^^^ ■ Q^"^^), we may also expand P/(n • q^"^^) in the expression for /^^^ in Eq. (35), 
using the addition theorem of spherical harmonics: 

m=—l 



Substituting Eqs. (35) and (37) into Eq. (34), and using the orthogonality of spherical 
harmonics, we obtain 

j> d'nf'\n ■ ^-W) X /(^)(n • g (^)) = £(2/ + 1) ttPK?^^^ • Q^'^) f['' fP- (38) 



Now, if either f'^^^ or f^'^^ is approximately band-limited to Lcgcndre polynomials with 
/ < /q, then the sum in Eq. (38) is over a finite number of terms. We shall assume the 
approximate band-limitcdncss of the jet profile e{n-b) in what follows. For jets characterized 
by an opening half- angle 9q, the band limit lo will be in excess of tt/Oq. 

The usefulness of Eq. (38) hinges upon our ability to calculate the Legendre transform 
integrals of Eq. (36), and to do so accurately and efficiently. We now consider how this is to 
be accomplished. 

It is in principle possible to use the well-known recursion relations for the Legendre 
polynomials to obtain expressions for the higher-/ transforms in terms of lower-/ transforms, 
at least for the transforms of the Doppler factor and for the Fisher jet profile. Such a 
recursion scheme is unquestionably efficient. Unfortunately, it turns out to be inaccurate — 
the resulting recursion relations are unstable, and for physically interesting values of /q the 
transforms at the largest values of / can be dominated by amplified numerical noise. 

The "Fast Lcgcndre Transform" (FLT) algorithm of Driscoll & Hcaly (1994) is well- 
suited to our purposes. It is the spherical analogue of the celebrated FFT of Cooley & 
Tukey (1965), and has already become somewhat well-known in astrophysical circles, among 
analysts of CMB data (Oh, Spergel, & Hinshaw 1998; Wandelt & Gorski 2001) 

Driscoll & Healy (1994) showed that for band-limited functions /*-*-'(x) with band limit 
lo, the Legendre transform of Eq. (36) is exactly given by the quadrature formula 

2^0-1 

/»= ^4'»V^^Hcos^,)PKcos^,), (39) 

j=0 



where 
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In Eq. (41) we have corrected the normahzation in the expression for the weights of DriscoU 
& Healy (1994), which was missing a factor \/2. 

Eq. (39) is the Discrete Legendre Transform (DLT), which is entirely analogous to the 
Discrete Fourier Transform. The band-limit /q plays a role analogous to that the Nyquist 
frequency - any power in the Legendre spectrum of f^'^\x) above Iq is aliased by the DLT 
down to harmonics below Iq. For this reason, computation by means of spherical expansion 
is suitable for smooth jet profiles, such as the Fisher profile or the power-law profile. Profiles 
with sharp edges, such as the top-hat jet profile, can be expected to have substantial power 
at very high harmonics, and are therefore less well-suited to computation by this method. 
We will exhibit compact analytic emission formulas for the top-hat jet profile in §3.3 below. 

Naive evaluation of the DLT evidently results in an expensive algorithm, one that is 
O{lo^). Physically interesting GRB jets are believed to be characterized by beaming angles 

~ 10~^ — 10~^ radians, and by jet opening angles > 10~^ radians, implying values of 
Iq in the several hundreds, or even above 1000. Particularly if many integrals are to be 
evaluated (as in a population-synthesis simulation, for example), the naive algorithm exacts 
a prohibitive cost. 

Fortunately, DriscoU & Healy (1994) also exhibited the FLT: a fast — 0(^0 log(Zo)^) — 
algorithm modeled after the "Divide and Conquer" strategy of the FFT (Cooley Sz Tukey 
1965). We have implemented this algorithm as a C library, and find it quite satisfactory for 
our present purposes. All computations of emission from Fisher and power-law profile jets 
exhibited in the next section were performed using this method. 



In §3.1 we saw that the spherical convolutions in which we are interested combine e{n-h) 
or (n ■ h)e{n ■ b) with (1 - /3n • f)"" {n = 2, 3, 4), or with {n ■ V){1 - f3n ■ 0"^ A DLT (that 
is, a list of Iq numbers) is required for each of the above convolvees. If each such DLT were 
computed as a separate FLT, the computational burden would rise to near-irritating levels. 

Fortunately, it is not necessary to compute an FLT for each of these DLTs. The well- 



3.2.2. Savings 
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known recursion relation {21 + l)xPi{x) — {I + l)Pi+i{x) + lPi-i{x) among the Legendre 
Polynomials, together with the relation (1 — f3x)~'^~^^ = (1 — f3x) x (1 — f3x)~'^ allow us to 
relate the DLT of (1 -/^x)"" (n = 2,3), and the DLT of a;(l-/3x)-^ to the DLT of (l-/?x)-^ 
and also to relate the DLT of x x e{x) to the DLT of e{x). 

It is thus only necessary, for the purposes of the present work, to calculate the DLTs of 

e{x) and of (1 — /5a;) The other DLTs are then computable essentially for free. If it were 
not of interest to calculate the distributions of §3.1.3, it would only be necessary to compute 
the DLTs of e{x) and of (1 — Px)~^, in order to compute Eis^ and ^'peak- 

Note also that since the two required DLTs are independent of the viewing angle 9^, the 
processing expense of calculating those DLTs is incurred once, after which the stored DLTs 
may be reused for many values of the viewing angle 9^. 

3.2.3. Numerical Limitations 

Eq. (38) has the form of an inner product with a weighting function proportional to 

Pi{q^^'^ ■ g (2)). In the cases we consider here, the argument of the Legendre polynomial is 
the cosine of the viewing angle, cos^^. When 9,^ is close to zero, the Legendre polynomial is 
nearly 1, and the terms in the inner product add constructively. This behavior is necessary 
so that maximum emission is obtained when the observing direction is aligned with the jet 
axis. 

As the observing direction moves away from the jet axis, terms in the inner product 
must begin to engage in some mutual cancellation, in order to produce the expected decrease 
in fiuence. If the underlying profile spans many orders of magnitude (as the Fisher profile 
may easily do), then it may come to pass that the cancellations required to produce the 
necessary suppression of fiuence at large 9^ cannot occur in practical computations, as they 
are drowned out by numerical noise due to the finite machine precision. 

The most commonly encountered form of fioating-point double precision is constructed 
from 8-byte words, which afford a relative spacing between adjacent representable floating 
point numbers of about 10~^^. This is the smallest relative magnitude that can be attained, 
using such arithmetic, in calculations that depend on delicate cancellation. 

If we demand a minimum accuracy of 1%, we may therefore expect that the computable 
dynamic range of Ei^o (say) is no more than about 12 decades, using 8-byte reals. This 
expectation is in fact borne out by experience. Of course this is plenty of dynamic range for 
the present purposes. Should greater dynamic range be nonetheless required, some additional 
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decades may be secured by adopting extended precision (typically 16-byte floating point 
words), and possibly computing with machine architectures endowed with a larger natural 
word size. 



3.2.4- "Smooth" Power-Law Profiles 

Power-law profiles diverge at the jet axis. It is customary to deal with this deplorable 
behavior by truncating the power- law at some small angle (Zhang & Meszaros 2002). This 
practice represents a problem for our spherical-expansion approach, since the sharp edge at 
the cutoff injects unwanted spectral power into the DLT of epL at high values of I. 

Consequently, the normalized power-law profile that we actually use in this work is 

epL(cos 9) = n{l -cos9 + Oq^ (42) 

where we have introduced a smooth cut-off ^o, rather than truncate the power- law. The 
resulting profile has a smooth, cusp-free peak in ^ at ^ = 0. 

The (single-lobe) normalization factor k is given by 

2^[(2+eoV2)l-^-(eoV2)l-^J ■ ^ (43) 
27rlog(46»o"^+l) ■ ^ 

The effective angular width of this profile is a function of both and 5. A simple 
estimator of this width may be constructed by making a Gaussian approximation to the jet 
near the peak. The width of this Gaussian is then approximately 



9yj = 



= eo/V2S. (44) 

Thus for a reasonable range of d near S — 1, 9q is itself not a bad estimator of the jet width. 
For large S, however, it is advisable to regard 9^;, rather than 9o, as the effective angular 
width of the profile. 

In particular, when we calculate E^^ = Eiso x (1 — cos9br) for power-law profiles, 
we always use a break angle 9br — ma'x.{9y,9w), rather than the "naive" prescription 9 — 
msix{9y, 9q). 
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3.2.5. Emission From Afterglow Jets 



The main focus of the current work is prompt GRB emission, so we are not really 
concerned with afterglow emission. Nonetheless, it is worthwhile to digress briefly on the 
applicability of our spherical convolution strategy to afterglows. 

In prompt emission, it is quite reasonable to assume that 7 is constant throughout the 
duration of the burst (Rees & Meszaros 1994). This assumption underlies our expressions 
for the various emission formulas derived above. 

Afterglow emission is another story. The external shock model features a deceleration 
of the shock front, due to interaction with the ISM (Piran 1999). This means that photons 
received simultaneously by the observer from different parts of the shock front were actually 
emitted at epochs characterized by different values of 7 — there is a relativistic "look-back" 
effect that must be taken properly into account (Woods & Loeb 1999; Granot et al. 1999). 
It may seem that this effect compromises the usefulness of the methods introduced in this 
section, by making the main emission formulas take forms that cannot be cast as convolutions 
on the 2-sphere. 

In fact, subject to the assumption that the shape of the underlying emission profile 
is invariant (so that only the emission amplitude varies with time), it is still possible to 
write the emission formulas as 2-sphere convolutions, so that the present approach using 
decompositions in Legendre polynomials may still be valuable for afterglow studies. 

This is the case because at each observer time tg, a photon received from the element of 

— * — * — * 

the jet at n was emitted at a time t{to, n-l) that is a function of n • Z, where I is the viewing 
direction. This function is computable in terms of the jet velocity history, /3(i). If the jet 
emission profile (now a function of time) has by assumption the form 



then the instantaneous flux is proportional to a convolution like that of Eq. (34), where 



It follows that under a moderately restrictive assumption about the time development of 
the emission profile, afterglow emission may also be treated by decomposition into Legendre 
polynomials, as described above. A possible generalization of this approach that relaxes 
the above restriction somewhat is to consider profiles constituted of linear superpositions of 
terms like the one in Eq. (45). 




(45) 




(46) 
(47) 
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Note another assumption made here, however: we still require that 7 is uniform across 
the jet. In fact, afterglow jet models with spatial variation in 7 are considered in the literature 
(Rossi ct al. 2002; Granot & Kumar 2003). If the uniform-7 assumption is removed, the 
required factorization of the integrand into and f^'^\n-l) fails, and direct expansion 

in Legendre polynomials is no longer an available strategy. If the variation of 7 across the 
jet is not too strong, however, it is still possible to adopt the spherical expansion method 
after expanding the Doppler factor in powers of the deviation of 7 from its value on the jet 
axis, 7(71) — 7(&). 



3.3. Emission From Top-Hat Jet Profiles 

As mentioned above, the sharp edge of a top-hat jet profile forces unwanted power onto 
the higher harmonics of its DLT. Fortunately, there are closed-form analytical expressions 
for the emission from a top-hat jet profile, which we present here. 



3.3.1. Emission Formulas 

The emission profile of a top-hat jet of opening half- angle 9q is 

eTHicosO) = [27r(l - cos ^0)]"^ ©(cos ^ - cos^o), (48) 
where Q{x) is the Heaviside step function. Prom Eq. (23) we then have 



27r74(l -cos^o) / (l-/3n-/)3 



E?^ 



dcos9 / — 3, (49) 



2717^(1 - cos 60) Jcoseo J-n [1 - /9(cos 6y COS 6 + sin 6^ sin 6 cos 0)] 

— * — * 

where again 9^ is the angle between b and I. 

The integral over (p may be performed using formula 3.661.4 of Gradshteyn & Rhyzik 
(1965, p. 383). The remaining integral over cos^ may then be performed by means of the 
successive substitutions y — (5 cos 9 — cos^„, tan 5 — lul sin^„. The result is 

2p7^(l — cosb'o) 

where 

u ^ _ 7^(27' - + (37' sin' 9, -l)z + 2 cos9, sin^ 9, 

(^2 + ^-2sin2e,)3/2 ■ 



-22- 



Prom Eq. (26), the average Doppler factor (D) is expressed as a ratio of two integrals, 
one of which is the same as the one in Eq.(49), while the other is quite similar and may be 
calculated using the same technique. The result is 

^ _i fiP - cos 0y) - f{(5 cos Oq - cos e^) 
^ ' ^ g{(5 - cos 9^) -g{(5 cos Oq- cos 9^)' ^ ' 

where f\z) is as in Eq. (51), while the function g{z) is given by 

272^ + 2 cos 
(2;^ + 7 ^ sm 9y)'^'^ 

The frequency distribution of Eiso is generally given by Eq. (27). This formula may be 
applied with Eq. (50) to yield an expression for dN/dEiso that is suitable for computation. 
This expression is algebraically obvious, but burdensome, and will be omitted here. 



4. Jet Profiles And Their Signatures 

We now turn directly to the observable signatures of universal jet models based on 
Fisher, power-law, and top-hat jet profiles. 

Our approach is to explore the parameter space of universal jet models based on these 
families of profiles, assessing features of the landscape for their resemblance or dissimilarity to 
the empirically-known properties of GRBs discussed earlier. We do not intend here to make 
quantitative judgments of how well the models "fit" the data. Rather, we wish to construct 
a map of the parameter space of universal-jet GRB models, so as to clarify the strengths 
and weaknesses of each profile choice. This allows one to choose which aspects of the data 
can be modeled successfully using a particular family of universal-jet models, and which 
must be "evaded" by invoking new model features. It also allows easier interpretation of the 
results of population-synthesis simulations that attempt to reproduce the distributional and 
correlational properties of E'iso, E^^\ and E^^ak, irrespective of whether the models adopted 
are of the universal or variable-geometry variety. 

This map is also a useful tool with respect to models that make more realistic observable 
choices than ours. Our variables are admittedly simplistic, insofar as comparison with data 
from real GRB detectors are concerned, i^iso is a bolometric quantity, which makes no 
reference to detector bandpass, and {D) x e'^^^ is a somewhat crude representation of the 
true spectral peak. More realistic modeling ascribes an actual rest-frame spectrum to the 
jet, possibly even allowing for spectral variation over the surface of the shock. Even more 
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realistic would be the introduction of time variability, so as to allow more detailed modeling 
of detector thresholds. 

What such modeling gains in fidelity, however, it loses in performance. It is an extremely 
burdensome task to compute emission, correlation, and distributional properties from models 
that must integrate numerically across the shock surface, the photon frequency domain, and 
the time domain. Such investigations cannot currently be carried out over wide swaths of 
parameter space. This is where our map can be valuable, by indicating qualitatively the 
kinds of models and relatively narrow ranges of parameters that are worth a closer look, and 
where builders of complicated models might consider unlimbering their heavy machinery. 

In this section we consider the following families of universal GRB jet models: 

• Fisher profile universal jets 

• Power-law profile universal jets, with power-law indices S—1, 2, and 8. 

• Top-hat profile universal jets. 

For each family of models, we exhibit six kinds of plots: Eiso versus 6^, (D) versus 6.^, 
E\^a versus {D) ("Amati" plot), versus {D) ("Ghirlanda" plot), dN/d log Eiso frequency 
plot, and dN/dlogE^^ frequency plot. 

In each case we allow 7 to take the values 10, 33.3, 100, and 333. 

Each family of profiles has a parameter which characterizes the angular width of the 
profile. Its interpretation is specific to the functional form of the profile. In almost all cases, 
we allow 9o to take the values 0.1, 0.3, 0.01, and 0.003. The one exception is the power-law 
profile with 6 = 8. In this case we actually select Oq = 0.4, 0.12, 0.04, and 0.012, so as to 
make the effective width estimator 9^ of Eq. (44) take the values 9yj =0.1, 0.3, 0.01, and 
0.003. 

The region of parameter space covered by these choices of parameter values includes 
a subregion where 7"^ is larger than the characteristic width of the underlying jet profile 
(^0 oi' ^wj) the case may be). This subregion of parameter space may be considered 
unphysical on hydrodynamic grounds (Rhoads 1999; Sari et al. 1999), since a hydrodynamic 
jet expands sideways rather than expanding ballisticaliy with a constant shape in this regime. 
Furthermore, general arguments of hydrodynamic causality during the acceleration phase 
suggest that the angular scale of the jet should be quite comparable to 7"^ (Stern 2003, 
Lazzati 2005, private communication). On the other hand, there is evidence that MHD jets 
can sustain degrees of coUimation that easily violate this constraint (Lapenta & Kronberg 



2005). We do not intend to discuss the details of these physical pictures here. We simply 
consider this region in our study for the sake of completeness. 



4.1. 



Emission Prom Fisher Jet Profiles 



In what follows, we use the single-lobe normahzed Fisher jet profile: 



eFisher(cOs6l) 



(54) 



'0 



4.1.1 



Emission 



The four panels of Fig. 3 show the behavior of Eiso as a function of viewing angle 6^. The 
noteworthy feature of these plots is that the value of 7^^ sets a lower limit to the effective 
angular size of the emission pattern. The characteristic angular width of radiation patterns 
from jets with 9o < 7^^ is set by 7^^, not 6*0, since in this limit the jet effectively resembles 
"needle jet," whose emission pattern is essentially set by the Doppler function. This effect is 
particularly conspicuous in the lower-right hand panel of Fig. 3, in which the intrinsic profile 
width is 0.003 radians, while all of the emission patterns are in fact considerably broader 
than this. 

This jet-broadening effect is a ubiquitous feature of relativistic kinematics, that recurs 
in all the other profiles as well. 

The four panels of Fig. 4 show the behavior of (D) as a function of viewing angle 
9y. Here a curious feature requires explanation: there is a "shelf" that appears in the (D) 
function near the jet axis, which for large 7 extends considerably beyond the value of ^o- 
For ^0 = 0.1, the 7 = 333 curve has a shelf that extends out to 6^ — 0.5, taking up about 
12% of a hemisphere. 

The explanation of the shelf is seen in the four panels of Fig. 5. These plots show 
the photon "provenance" distribution for an observer viewing a Fisher jet with parameters 
7 = 333, ^0 — 0.1, at viewing angle 9^. The curve in each panel is a slice through the plane 
containing I and b of the function eFishcr(^ ■ ^) x (1 ~ /^^ ■ ^)^^! plotted as a function of angle 9 

— * 

from b. Ignoring the azimuthal direction for the sake of simplicity, this is the (un-normalized) 
probability that a photon directed along / should have been emitted from a location on the 
shock surface that is at an angle 9 from the jet axis. The four panels show this curve for four 
different viewing angles — 9y = 0.2, 0.45, 0.6, and 0.75. The corresponding point on the 
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corresponding (D) — 9^ plot of Fig. 4 is shown in the inset. Angular scales are represented 
linearly in these plots, the better to bring out the detail. 

The dashed line shows the position of the "Doppler-mean" viewing angle {9) — that is 
the angle such that 7~^[1 — /3cos{9v — {9))]~^ — (D). It represents the provenance of the 
photons most responsible for giving (D) its current value. 

It can be seen from these plots that (^^) is simply tracking the mass of the distribution, 
whose location and extent is determined by the interplay between the Fisher function and 
the Doppler photon distribution. Near the axis of the jet, the peak of the Doppler factor 
dominates, and almost all the received photons originate from locations on the jet that are 
moving towards the observer (top left panel). 

When the viewing angle gets large enough, the Doppler peak starts to peel away from 
the jet axis. The provenance distribution begins to develop a bimodal character (with the 
Doppler mode still dominant), so the mean photon provenance shifts slightly away from the 
part of the jet moving towards the observer, towards the axis of the jet. This is the beginning 
of the shelf drop-off (top right panel) . 

As the viewing angle continues to grow, the jet axis peak of the distribution becomes 
more and more competitive with the Doppler peak, pulling (9) substantially away from the 
part of the jet moving towards the observer (bottom left panel). 

Finally, at large viewing angles, the Doppler peak is eclipsed by the jet axis peak. The 
value of (9) approaches ^o, and changes very little thereafter. From here on out, the jet 
appears to the observer as a needle, with no angular structure as far as {D) is concerned 
(bottom right panel). 

The presence of this shelf at the bright end of universal Fisher jet models with 7""^ <S ^0 
is potentially very constraining on such models. We have already noted that in the case of 
7 = 333, 6*0 = 0.1, the shelf takes up 12% of a hemisphere. In the context of a universal 
model, this means that for these parameters, the brightest 12% of GRBs should belong to 
a sub-population with constant (D) and hence with constant -Epcak- Such a concentration 
might very well be conspicuous. As we will see, it also strongly distorts the £'iso-(-D) and 
Ej^-{D) relations away from their empirical forms. 

4.1.2. Eiso - {D) and E^^^ - (D) Correlations 

The four panels of Fig. 6 show the behavior of the E[so^{D) ("Amati") correlation in 
universal Fisher models with our standard parameter choices. The axis scales are chosen 
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so that a slope of 0.5 appears as a 45° line, in order to ease the search for an "Amati-like" 
slope. 

The 1/3 slope characteristic of off-axis emission from "needle" jets is seen to hold for 
^0 <^ 7~^- This part of parameter space is unpromising, since the empirical slope is known 
to be 0.5 (Amati et al. 2002; Lamb et al. 2004). 

For <^ ^0, there does appear a section of the Eiso-{D) curve that seems to have the 
empirically-observed slope. This part of the curve does not, however, include the brightest — 
and most easily observable — events. In fact, the "Amati-like" part of the curve corresponds 
to the narrow "cliff" of Fig. 4. For 7 = 333, 9o — 0.1 (lower left panel), for example, it 
comprises only viewing angles 9y between 0.5 and 0.7, subtending only 11% of a hemisphere, 
compared to the 12% subtended by the shelf. The case of 7 = 333, 9q = 0.03 (upper right 
panel) seems more promising, since now the shelf sets in at about 9y = 0.1, subtending a 
much smaller solid angle. 

It appears, then, that a fit of a universal Fisher jet to "Amati"-type data could be 
extremely constraining on the model parameters. The necessity of producing the Amati 
slope, while at the same time not inflicting on the model predictions a sub-population of 
events from the shelf with zero slope that is brighter than, and comparably numerous to, the 
"Amati-like" sub-population, would in all likelihood result in a narrow acceptable region of 
7-^0 parameter space. 

The four panels of Fig. 7 show the behavior of the E^^^-{D) ( "Ghirlanda" ) correlation 
in universal Fisher models with our standard parameter choices. 

The < limit appears to feature a slope of 0.5, as opposed to the empirically- 
determined value of 0.7 (Ghirlanda et al. 2004), so again this part of parameter space is not 
promising territory for this family of models. 

For 7~^ -C ^0; there is a steeper slope, corresponding again to the drop-off from the 
shelf (as may be determined by comparison with Fig. 4). For = 0.1, 7 = 333 (upper left 
panel), this slope is 1. It drops down to zero at high values of E'^\ Similar behavior is 
seen for 9q = 0.03, 7 = 333 (upper right panel). The difference between the two curves is 
that the flattening-out takes place over a narrower range in 9^ when ^0 = 0.03 than when 
9q — 0.1, reducing the anomalous zero-slope sub-population. So just as for the Amati case, 
it appears that a fit of a universal Fisher jet to "Ghirlanda" -type data could produce very 
strict constraints on 7 and ^o- 

One feature of Fig. 7 requiring an explanation is the hook-like structure observed for 
7~^ > ^0, in which e!^^ appears to fold back upon itself. 
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This feature is another side-effect of the jet-broadening effect noted above. The rela- 
tivistic kinematics is cutting off the jet profile when 6^ ~ 7"^, long before 6^ approaches ^o- 
Therefore, the factor 1 — /it, in the definition of E^'^ (Eq. 10) begins to depress the value of 
e'^'' earlier and more deeply than it could the "bare" fisher profile. With reference to the 
schematic Fig. 2, the maximum of the curve is occurring at lower values of fiy, and the curve 
plunges to much lower values before it " recovers" at 0^ — ^o- 

What this indicates is that when relativistic kinematics are taken into account, it would 
be best to apply the Frail et al. (2001) procedure using a break angle 6*6^ = max(^^, ^0, 7~^)) 
in order to avoid such artifacts. 



4.1.3. Eiso and Ey Frequency Distributions 

The four panels of Fig. 8 show the effect of relativistic kinematics on the uniform E[so 
frequency distributions of universal Fisher jets that were derived in §2.1. In each panel, the 
light solid line shows the "bare" value of the distribution, corresponding to the (uniform) 
shape derived just from the profile. Essentially, this is the 7 — s> 00 limit. 

We see here that with 7"^ -C 9o, the uniform distribution that was built into the Fisher 
profile by construction is largely preserved for bright events. As 6q decreases, the jet emission 
behaves more and more like off-axis emission from a needle jet, which by Eq. (32) is a power- 
law with a slope of -1/3. Accordingly, the frequency distributions depart more and more 
from the uniform shape that was due to the Fisher profile. 

The four panels of Fig. 9 show the effect of relativistic kinematics on the broad E^ 
distributions of universal Fisher jets that were derived in §2.1. The "bare" distribution is 
again shown by the light solid line. Just as for the ii^iso frequency plots, these distributions 
remain broad, in resemblance to the "bare" version, until ^0 becomes so small that the 
emission is characteristic of an off- axis needle jet, at which point they assume the shape of 
a power-law with a slope of -1/2, as expected from Eq. (33). 



4.2. 5 — 1 Power-Law Jet Profiles 

4.2.1. Emission 

The four panels of Fig. 10 show the behavior of E'iso as a function of viewing angle 9y for 
the power-law profile with index S = 1. Once again, we sec the "jet-broadening" effect, in 
which the effective angular size of the jet radiation pattern is set by whichever is larger of ^0 
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and 7 ^. These radiation patterns have broad tails in every case, reflecting the fact that this 
profile falls off more gradually than the Doppler photon distribution function 7~^(1 — 

The four panels of Fig. 11 show the behavior of (D) as a function of viewing angle 9y. 
The striking feature of these plots is the extremely restricted dynamic range of (D) . 

The large-6'^ limit in each case is (Z^) = 7, as is expected by the fact that far from 
the jet axis, this profile falls so gradually that it looks nearly uniform (recall from §3.1.2 
that 7 is the "uniform outflow" limit of (D)). This reason for this large-^^ hmit of (D) also 
helps explain why the dynamic range of (D) is so limited — this profile varies quite gently 
almost everywhere, so its (D) deviates only moderately from 7 almost everywhere. The only 
real chance that (D) has to run away is near 9^ = 0, but even here it is eventually reined 
in by the cutoff. The smaller 6q, the wider the available range of (D), but it is clear that 
uncomfortably small values of 9q will be required in order to allow this profile to produce 
the kind of range of £'peak — at least two orders of magnitude — that is observed in GRBs. 

4.2.2. Ei,o - (D) and E^^^^ - (D) Correlations 

The four panels of Fig. 12 show the behavior of the Eiso^{D) ("Amati") correlation in 
universal S = 1 power-law models with our standard parameter choices. 

Here we see the some of the trouble created by the limited dynamical range of (D) . The 
curves stay more-or-less flat for much of the dynamical range of Ejso- The high-7 curves do 
suddenly rear up at high i^iso, but not nearly enough to make a convincing Amati plot, and 
in any event, comparison with Figs. 11 and 10 show that this is going on over a very small 
region of 9^, corresponding to a tiny fraction of observed bursts in a universal jet model. 

The four panels of Fig. 13 show the behavior of the Ef^^-{D) ( "Ghirlanda" ) correlation 
in universal S — 1 power-law models with our standard parameter choices. 

There seems to be no hope for a real "Ghirlanda" relation with this profile. Due to the 
limited dynamic range of (D), the trajectory is essentially fiat, with a hook-like structure 
at the very end due to the same side-effect of relativistic jet-broadening that was remarked 
upon in §4.1.2. 

4.2.3. Eiso and E^l^ Frequency Distributions 

The four panels of Fig. 14 show the effect of relativistic kinematics on the E^so frequency 
distribution ol5 — 1 power-law jets — see Eq. (18). Again, the "bare" (7 — > 00) distribution 
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is shown by a light sohd hne, whose slope in this case is -1. 

We see here that the distribution of dN/dEiso appears to be irreducibly narrow, and 
relativistic kinematics are powerless to broaden it. Prom this it appears that this "extremal" 
power-law profile cannot under any circumstances reproduce the broad — and very nearly 
uniform in log E'iso — distribution that is observed empirically in GRBs. 

The four panels of Fig. 15 show the relativistic broadening of the 5-function distribution 
designed into this profile — see Eq. (19). It clearly succeeds in preserving its empirically- 
desirable narrowness. 

4.3. 5 — 2 Power-Law Jet Profiles 

The results of §4.2 are to some extent anomalous for power-law jets, since 5 = 1 is a 
singular case as far as Ej^ is concerned. The 5 — 2 case comprises more typical behavior. 

4-3.1. Emission 

The four panels of Fig. 16 show the behavior of Eiso as a function of viewing angle 6^ for 
a power-law profile with 6 = 2. Wc sec jet broadening again, as expected. The tails at high 
6y are no longer independent of 7, as they were for 5=1, since now the profile function falls 
off as rapidly as the Doppler photon distribution function. The available dynamic range of 
Eiso is thereby increased. 

The four panels of Fig. 17 show the behavior of (D) as a function of viewing angle 6y. 
We see that (D) can now just about produce two orders of magnitude's worth of dynamic 
range, in the "needle jet" approximation (^0 <S 7~^)) but the dynamic range again becomes 
quite restricted in other regimes. For 9o > 7~^, the large-^^ limit of (D) is still 7, as expected, 
and this hmit still serves to constrain the available range of (D). 

4.3.2. Eiso - {D) and E^f^ - {D) CorrelaUons 

The four panels of Fig. 18 show the behavior of the Eiso-{D) ("Amati") correlation in 
universal 5 = 2 power-law models with our standard parameter choices. 

Here it appears that this profile is trapped between a rock and a hard place. In order to 
make enough dynamic range in (D), it must enter the "needle jet" regime (^0 7~^)- But 
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as soon as it does so, the slope of the Eiso-{D) curve becomes 1/3, as expected for off-axis 
emission from narrow jets. It does not appear that this profile can make an "Amati" relation 
spanning two decades of Ep^ak- 

The four panels of Fig. 19 show the behavior of the E^^-{D) ( "Ghirlanda" ) correlation 
in universal 5 — 2 power-law models with our standard parameter choices. 

We can see from these curves that this profile's abihty to make a respectable "Ghirlanda" 
relation is no better than it's ability to make "Amati" relations, and for the same reasons. 

4. 3. 3. Eiso o-nd E^^^ Frequency Distributions 

The four panels of Fig. 20 show the effect of rclativistic kinematics on the £^iso frequency 
distribution of 5 = 2 power- law jets — see Eq. (18). Again, the "bare" (7 — >• 00) distribution 
is shown by a light solid line, whose slope in this case is -1/2. 

The distribution of dN / dEiso is still narrower than the data seem to indicate, although 
in the "needle-jet" limit 9q <C 7~^, the slope of the distribution is softened by rclativistic 
kinematics to the -1/3 index power-law expected from Eq. (32). 

The four panels of Fig. 21 show the effect of rclativistic kinematics on the Ej'^ frequency 
distribution of power-law jets — see Eq. (19). The distribution can easily be made acceptably 
narrow, although in the "needle-jet" limit 6*0 -C 7"^, the -1/2 slope expected from Eq. (33) 
seems somewhat less steep than is consistent with the data. This means that a fit of such a 
universal model could potentially provide a (model-dependent) lower limit on 7. 

4.4. 5 — 8 Power-Law Jet Profiles 

The chief difficulty encountered by the relatively shallow 6 = 2 index power-law of §4.3 is 
a lack of dynamic range, particularly for (D). It is therefore interesting to examine whether 
this weakness can be addressed by adopting a steeper power-law index. Here we examine 
5 = 8. As a reminder, our standard parameter values are altered here, to ^0 — 0.4, 0.12, 
0.04, and 0.012, so as to make the effective width estimator 9^ of Eq. (44) take the values 
Oyj =0.1, 0.3, 0.01, and 0.003. 
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4.4- 1- Emission 

The four panels of Fig. 22 show the behavior of i^iso as a function of viewing angle 9^ for 
a power-law profile with 6 = 8. We see jet broadening again, as expected. For larger values 
of 7, this profile appears to be able to produce a fairly wide variety of radiation patterns, 
and great dynamic range of £^iso- 

The four panels of Fig. 23 show the behavior of (D) as a function of viewing angle 
6y. Here we can finally see strong dynamic range of {D) in a power- law model, at least for 
narrow jets. We also see the return of the "shelf" feature that was noted in §4.1.1, which 
appears to be a generic feature of steep profiles. 



4-4-2- Eiso - (D) and Ey - (D) Correlations 

The four panels of Fig. 24 show the behavior of the Eiso-{D) ("Amati") correlation in 
universal 6 = 8 power-law models. The axis scales are chosen so that a slope of 0.5 appears 
as a 45° line, in order to ease the search for an "Amati-like" slope. 

Here it appears that slopes steeper than 1/3 can be produced, possibly even over a 
sufficiently large dynamic range of i^iso and (D), at the cost of some fiattening out near 
the bright end of the distribution. The slope of the 6o = 0.04, 7"^ = 0.003 curve (bottom 
left panel), for example, is about 0.68 at intermediate ii^iso, and by comparison with Fig. 23 
(bottom right panel), the flattening out at high values of (D) only sets in at about = 0.1, 
subtending only 5 x 10~^ of a hemisphere. While the slope is high compared to the standard 
Amati value of 0.5, it appears possible that there could be a value of 5 bracketed between 
6 = 2 and 6 = 8 that can reproduce the Amati slope. 

The four panels of Fig. 25 show the behavior of the E^^-{D) ( "Ghirlanda" ) correlation 
in universal 6 = 8 power-law models. 

Just as for the "Amati" correlation plots, there appears to be a range of parameter 
space with relatively steep slope and large dynamic range of {D). The slope of the ^0 — 0-04, 
7-1 — 0.003 curve (bottom left panel), for example, is about 0.83 at intermediate Ei^o, and 
again the flattening at high energy subtends a tiny solid angle. It thus appears possible 
that there could be a value of 6 bracketed between 6 = 2 and 6 = 8 that can reproduce the 
Ghirlanda slope. 
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4.4-3. Eiso and Ey Frequency Distributions 

The four panels of Fig. 26 show the effect of relativistic kinematics on the £^iso frequency 
distribution of 5 = 8 power-law jets — see Eq. (18). Again, the "bare" (7 — > 00) distribution 
is shown by a light solid line, whose slope in this case is -1/8. 

The large value of S has finally made it possible for this distribution to become relatively 
shallow and broad, possibly doing less violence to the observed, nearly-uniform distribution 
than was possible using shallower power-law indexes. It remains the case that the "needle 
jet" regime 6*0 <^ must be avoided to avoid introducing too great a non-uniformity (see 
Eq. (32). This is Hkely to place a lower bound on 7^0 in fits to the Eiso distribution data. 

The four panels of Fig. 27 show the effect of relativistic kinematics on the e!^^ frequency 
distribution of 5 = 8 power-law jets — Eq. (19). Again, the "bare" (7 — > 00) distribution is 
shown by a light solid Hue, whose slope in this case is -1/7. 

Here the case is reversed. The "bare" slope is far too shallow, compared to the strongly- 
peaked distribution of the real data. It seems unlikely that relativistic kinematics in the 
"needle jet" regime can save the day, since even ignoring the just- discussed deleterious effect 
of this regime on the Eiso frequency distribution, the sharpest slope that can be produced 
in this regime is -1/2, by Eq. (33), which seems somewhat less steep than is consistent with 
the data. 



4.5. Emission Prom Top-Hat Jet Profiles 

4.5.1. Emission 

The four panels of Fig. 28 show the behavior of Eiso ^ a function of viewing angle 9^ 
for a top-hat profile. We see jet broadening again, as expected. 

The four panels of Fig. 29 show the behavior of (D) as a function of viewing angle 9^. 
There is no "shelf" effect here, in contrast to the Fisher and the 6 — S power-law cases — 
the sharp edge of the profile is the feature that dominates the emission for 7"^ < ^o- One 
curiosity is the small "horn" structure that appears at the edge of the profile when 7"-^ < ^o- 
For example, the curve for 60 = 0.03, 7^^ = 0.003 (top right panel) shows a brief upswing of 
(D) just before it plummets off the edge. This feature is easily explained: (D) represents an 
average of the Doppler factor weighted by the Doppler photon distribution function. Both 
the Doppler factor and the Doppler photon distribution function peak along the line-of-sight. 
Over most of the face of the profile, the value of the average Doppler factor is dragged down 
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from the value at the peak by the lower values corresponding to photons arriving from off the 
line-of-sight. However, within 7"^ of the edge of the jet, some of the photons that perform 
this "dragging down" office are simply missing, so the average briefly surges up. 

4.5.2. Eiso - {D) and Ej^ - {D) Correlations 

The four panels of Fig. 30 show the behavior of the Eiso-{D) ( "Amati" ) correlation in 
universal top-hat models. The axis scales are chosen so that a slope of 0.5 appears as a 45° 
hne, in order to ease the search for an "Amati-hke" slope. 

As expected, the "needle jet" cases 9o <S 7"^ exhibit the 1/3 slope characteristic of 
off-axis emission. These cases are therefore unpromising. The opposite case, 9o > 7~^ is 
more interesting. For example, in the case ^0 = 0.1, 7"^ = 0.003 (top left panel), the slope 

of the curve near the top is 0.46, and the gentler slope only sets in for fainter events. The 
similarity of the case 6q = 0.03, 7 = 333 (top right panel) to Figure 1 of Toma et al. (2005) 
(which assumed Oq = 0.02, 7^^ = 300 is quite striking, particularly in view of of the contrast 
between the rather detailed nature of the Toma et al. (2005) model and the rather simple 
model presented here. 

Potentially troubling is the fact that all the events seen on-axis are located at the tip 
of the hook structure at the end of the curve (this structure is a side-effect of the "horn" 
structure remarked upon in the previous section). Therefore, fits of this model to "Amati" - 
type data must constrain 9q to be small enough as not to overload the top end of the Amati 
relation with an excessive concentration of events. Note also that in order to take this profile 
seriously in a universal model, it is necessary to abandon the "GRBs are on-axis, XRFs are 
off-axis" paradigm of Yamazaki et al. (2002, 2003), since otherwise only XRFs would obey 
the Amati relation, while all GRBs would cluster very tightly in the £'iso--E'peak plane. 

The four panels of Fig. 31 show the behavior of the E^^-{D) ( "Ghirlanda" ) correlation 
in universal top-hat models. 

On the rising slope part of these curves, it appears that "Ghirlanda" -like slopes can 
indeed be made for 7"^ < ^o- As was the case for the Amati relation, the delicacy required 
to prevent a conspicuous blob of on-axis events from distorting the plot must necessarily 
constrain ^o- 
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4-5.3. Eiso and Ey Frequency Distributions 

The four panels of Fig. 32 show the effect of relativistic kinematics on the Eiso frequency 
distribution of top-hat jets. 

Here is evident the main difficulty that inheres in these models: the distribution of 
dN/ d log E[so simply can not be made even approximately uniform, as required by the data. 
Rather, for most of the dynamic range dN / dlogE-^^o ~ E-^^J , as expected for off-axis emis- 
sion (Eq. 32). The distribution then piles up at the top of the E^^^ range, where all the 
events that are viewed on-axis reside. It is apparent that no choice of parameters can make 
this model even remotely resemble the Ei^o frequency data. 

The four panels of Fig. 33 show the effect of relativistic kinematics on the e':^'' frequency 
distribution of top-hat jets. 

The curious, two-peak structure observed in some of these curves is produced by the 
application of the Frail et al. (2001) prescription to top-hat profiles, in conjunction with 
the discontinuous nature of the break angle, 9i„. = max{6^,6Q). A spike can appear near 
the energy corresponding to Oq, since the functional dependence of Ey'^ on 6^ can produce a 
maximum there, which produces a spike in the distribution function by Eq. (30). 

Events viewed near- or on-axis (the "spiky" components) appear capable of producing 
relatively narrow E^^ frequency distributions, especially for < ^o- The case of ^0 = 0.1, 

= 0.003 (top left panel) is an example. In this case, however, one must write off the 
off-axis component, which is undeniably broad, perhaps under cover of the fact that these 
would be XRFs, and cannot at present be placed in a "Ghirlanda"-type plot, since no break 
times have yet been observed for XRFs with redshifts. 



5. Discussion 

From the foregoing exploration, it appears that relativistic kinematics cannot alleviate 
the weakness of universal jet models that we identified in §2.4: it does not appear possible 
for a universal jet model to produce a narrow distribution of ^simultaneously with a broad 
distribution of E'iso- This is a difficulty for universal jets, although not necessarily a fatal 

one. The currently- available sample of E^"^ values is certainly more strongly affected by 
systematic selection effects than the sample of i^iso values, since there is some fraction of 
bursts for which no jet break is known. In particular, there are no jet breaks available for 
X-Ray Flash events, so that the frequency distribution of e!^^ may in fact extend to XRFs, 
despite being truncated by observational selection effects. Nonetheless, even concentrating 
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only on GRBs, it is striking that the empirical width of the distribution of GRBs is is 
not much larger than a single decade, while the width of the E'isodistribution is somewhere 
between three and four decades. It is hard to understand how such a discrepancy could be 
produced by selection effects alone. 

There is a generic effect of relativistic kinematics that is worth highhghting: The rel- 
ativistic beaming angle 7"^ imposes a lower limit on the effective angular size of the jet. 
Underlying jet emission profiles that are narrower than 7"^ will be broadened to an angular 
size of 7^^ by the relativistic Doppler effect. This effect has been noted in the case of after- 
glow emission (Woods & Loeb 1999; Rossi et al. 2002), but its influence on prompt emission 
has not been widely discussed. The net effect is to impose an upper bound on the extent to 
which the prompt flux can be amplified by coUimation of the profile. A concomitant effect 
is to curb the ability of very steep profiles to produce a large dynamic range of emission 
intensity. 

It seems worth commenting further upon the "shelf" feature of the (D) — 9y plots that 
was examined carefully in the case of Fisher jet profiles, and remarked upon in the case of 
5 — 8 power-law profiles. This feature appears to be a generic property of steep profiles, 
as might be expected from the explanation provided in §4.1.1. It illustrates the fact that 
relativistic kinematics can not only serve to soften a source spectrum (as is the case for 
the spectrum emitted by a narrow jet observed far off-axis), but can also harden a GRB 
spectrum with respect to what might be expected from a naive application of "needle-jet" 
kinematics. For example, in the case of a Fisher jet with 60 = 0.1, one might naively expect 
that for 9v > Oq, the jet profile presents the appearance of a needle jet, so that (D) must drop 
away as (1 — /3/i„)~^. It is therefore surprising to find that (D) can in fact remain relatively 
constant and considerably higher than this naive argument would lead us to believe, even 
quite far from the jet axis. 

It is encouraging to note that Amati-like and Ghirlanda-like correlations may emerge 
from "steep" profiles such as the Fisher profile, or the 5 — 8 power-law profile, under suitable 
choices of 9o and 7. It is somewhat striking, even surprising, to find such correlations 
emerging of their own accord from top-hat profiles. Such profiles are not usually considered 
promising candidates for universal models, as they are not "structured" . The main exception 
to this view has been put forward by Yamazaki and collaborators (Yamazaki et al. 2002, 2003, 
2004; Toma et al. 2005), as well as other papers by the same authors. It is certainly proven to 
be something of a surprise that top-hat models can produce Amati slopes much steeper than 
the 1/3 value expected from a relativistic point source (the "needle jet" approximation). 

It is noteworthy that in the case of the jet profiles that have a decent chance of making 
realistic-looking Amati/Ghirlanda relations (the Fisher profile, the 5 = 8 power- law profile. 
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and the top-hat profile), it is invariably the case that the Ghirlanda relation has a steeper 
slope than the Amati relation. This steepening is evidently a purely geometrical effect. It 
is clearly unrelated to the steepening effect noticed by Levinson & Eichler (2005), since we 
have certainly not included any dependency of the break angle on -Eiso, the source of the 
steepening noted in that paper. Apparently, with a break angle that is conceived purely 
in terms of geometry, it is still the case that Ghirlanda relations of universal models are 
expected to be steeper than their corresponding Amati relations. 

The inability of shallow power-law profiles (6 = 1 and 5 = 2) to produce any substan- 
tial dynamic range in (D) spells major trouble for such models, above and beyond their 
expected inability to produce broad Eisq distributions. This aspect of the situation might be 
correctable, if one is willing to incur the cost of introducing an angular dependence of the 
spectrum, in conjunction with the angular dependence of the emission profile. In particular, 
one could consider the (purely empirical) introduction of a relation like £^peak oc e{n- 6)^/^ as 
a microphysical relation satisfied on the surface of the shock. 

There have been some investigations that appear at first glance to do something like 
this. Yamazaki et al. (2004) introduced an intrinsic Amati relation between Ep^ak and E'lso 
in their universal off-axis model, to examine how the correlation was modified by relativistic 
kinematics for off-axis observers. Zhang et al. (2004) similarly assumed an intrinsic Amati 
relation in their study of Gaussian jet profiles. 

Kobayashi et al. (2002) pointed out that -Bpeak^ is proportional to the internal energy 
of the jet material. Zhang & Meszaros (2002b) and Ramirez-Ruiz & Lloyd-Ronning (2002) 
observed that internal synchrotron models can certainly produce a correlation between -Bpeak 
and the luminosity, although the exact shape of the correlation varies with the dependence 
of L on 7. It should be pointed out, however, that these dependencies of luminosity on 
-^pcak appear to be "global" statements about average conditions prevailing in the jet. The 
question of whether an internal shock model radiating by synchrotron emission might be 
expected to exhibit a local Amati- like correlation between Epeak and emissivity has not, to 
our knowledge, been addressed. 

We note in passing that it would be straightforward to modify our method to calculate 
E'peak instead of (D) in such a model, since all that is necessary is to introduce an additional 
multiplicative function Ep^akin ■ b) in the integral in the numerator of Eq. (26). 

We gratefully acknowledge Wayne Hu, Hiranya Peiris, Rosalba Pcrna, Davide Lazzati, 
and Ryo Yamazaki for detailed discussions and comments. This research was supported in 
part by NASA Contract NAGW-4690 and NASA Grant NAGS- 10759. 
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Fig. 1. — Empirical properties of E'iso and e!^\ Top left: Amati relation. Top right: 
Ghirlanda relation. Bottom left: Histogram of ^^iso values, excluding XRFs. Bottom right: 
Histogram of e!^'' values for events with well-determined The data used to make the 

plots in the left-hand panels is from Lamb et al. (2004). The data used to make the plots in 
the right-hand panels is from Ghirlanda et al. (2004). 
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essentially the 7 — > 00 limit. Its slope is -1, as expected from Eq. 18. 
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Fig. 15. — dN/d\og{Ej^) distributions for universal jet models based on S = 1 power- law 
profiles. Each panel shows the result for 7 = 10. 33.3, 100, and 333. Top left: = 0.1. Top 
right: 60 = 0.03. Bottom left: = 0.01. Bottom right: = 0.003. The vertical scale is 
arbitrary. The light solid line shows the "bare" distribution due to the underlying profile — 
essentially the 7 — > 00 limit. If it were not for the 9q cutoff and the relativistic kinematics, 
this distribution would be a 5-function, by construction. 
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Fig. 16. — Eiso as a function of viewing angle for emission from relativistic jets with 6 = 2 
power-law emission profiles. Each panel shows the result for 7"^ = 0.1, 0.03, 0.01, and 0.003. 
Top left: 60 = 0.1. Top right: = 0.03. Bottom left: 60 = 0.01. Bottom right: ^0 = 0.003. 
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Fig. 17. — Photon-number-averaged Doppler factor {D) as a function of viewing angle for 
emission from relativistic jets with 5 — 2 power-law emission profiles. Each panel shows the 
result for 7-^ = 0.1, 0.03, 0.01, and 0.003. Top left: = 0.1. Top right: Oq = 0.03. Bottom 
left: 9o = 0.01. Bottom right: = 0.003. The dynamic range of (D) is still quite limited, 
for almost all of the selected parameters. 
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Fig. 18. — Eiso-{D) ("Amati") plots for universal jet models based on 6 — 2 power- law- 
profiles. Each panel shows the result for 7-^ = 0.1, 0.03, 0.01, and 0.003. Top left: 60 = 0.1. 
Top right: ^0 = 0.03. Bottom left: ^0 = 0.01. Bottom right: ^0 = 0.003. 
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Fig. 19. — e'^^-IsJD) ( "Ghirlanda" ) plots for universal jet models based on 5 = 2 power-law 
profiles. Each panel shows the result for 7-^ = 0.1, 0.03, 0.01, and 0.003. Top left: = 0.1. 
Top right: ^0 = 0.03. Bottom left: = 0.01. Bottom right: ^0 = 0.003. 
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Fig. 20. — dN / dlog{Eiso) distributions for universal jet models based on 5 = 2 power-law 
profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: ^0 = 0.1. Top 
right: = 0.03. Bottom left: 60 = 0.01. Bottom right: = 0.003. The vertical scale is 
arbitrary. The light solid line shows the "bare" distribution due to the underlying profile — 
essentially the 7 — > 00 limit. Its slope is -1/2, as expected from Eq. 18. 
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Fig. 21. — dN / d\og{E'^^) distributions for universal jet models based on 5 = 2 power-law 
profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: 9q — 0.1. Top 
right: Oq = 0.03. Bottom left: Oq = 0.01. Bottom right: Oq = 0.003. The vertical scale is 
arbitrary. The hght sohd line shows the "bare" distribution due to the underlying profile — 
essentially the 7 — > 00 hmit. Its slope is -1, as expected from Eq. 19. 
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Fig. 22. — Eiso as a function of viewing angle for emission from relativistic jets with S — 8 
power-law emission profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top 
left: 00 = 0.4. Top right: 60 = 0.12. Bottom left: Oq = 0.04. Bottom right: Oo = 0.012. 
These values of are chosen so as to make the effective profile widths given by Eq. 44 be 
0.1, 0.03, 0.01, and 0.003, respectively. 
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Fig. 23. — Photon-number-averaged Doppler factor (D) as a function of viewing angle for 
emission from relativistic jets with 6 = 8 power-law emission profiles. Each panel shows the 
result for 7 = 10, 33.3, 100, and 333. Top left: ^0 = 0.4. Top right: 60 = 0.12. Bottom left: 
60 = 0.04. Bottom right: Oq = 0.012. 
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Fig. 24. — Eiso-{D) ("Amati") plots for universal jet models based on 5 — 8 power- law 
profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: 9o — 0.4. Top 
right: Oq = 0.12. Bottom left: 6*0 = 0.04. Bottom right: = 0.012. 
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Fig. 25. — e!^^-{D) ( "Ghirlanda" ) plots for universal jet models based on 5 — 2 power-law 
profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: 9o — 0.4. Top 
right: Oq = 0.12. Bottom left: ^0 = 0.04. Bottom right: ^0 = 0.012. 
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Fig. 26. — dN / dlog{Eiso) distributions for universal jet models based on 5 = 8 power- law- 
profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: Oq — 0.4. Top 
right: 9o = 0.12. Bottom left: Oq = 0.04. Bottom right: Oq = 0.012. The vertical scale is 
arbitrary. The light solid line shows the "bare" distribution due to the underlying profile — 
essentially the 7 ^ 00 limit. Its slope is -1/8, as expected from Eq. 18. 
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Fig. 27. — dN/d\og{E!f') distributions for universal jet models based on 5 = 8 power- law- 
profiles. Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: ^0 = 0.4. Top 
right: 9o = 0.12. Bottom left: 9o = 0.04. Bottom right: ^0 = 0.012. The vertical scale is 
arbitrary. The light solid line shows the "bare" distribution due to the underlying profile — 
essentially the 7 ^ 00 limit. Its slope is -1/7, as expected from Eq. 19. 
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Fig. 28. — £^iso as a function of viewing angle for emission from relativistic jets with top-hat 
emission profiles. These figures help illustrate the fact that 7"^ sets a lower hmit on the 
effective angular size of the jet. Each panel shows the result for 7"^ = 0.1, 0.03, 0.01, and 
0.003. Top left: 60 = 0.1. Top right: Oq = 0.03. Bottom left: Oq = 0.01. Bottom right: 
^0 = 0.003. 
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Fig. 29. — Photon-number-averaged Doppler factor {D) as a function of viewing angle for 
emission from rclativistic jets with top-hat emission profiles. Each panel shows the result 
for 7-1 = 0.1, 0.03, 0.01, and 0.003. Top left: Oq = 0.1. Top right: 60 = 0.03. Bottom left: 
00 = 0.01. Bottom right: = 0.003. 
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Fig. 30. — Eiso-{D) ( "Amati" ) plots for universal jet models based on top-hat profiles. Each 
panel shows the resuh for = 0.1, 0.03, 0.01, and 0.003. Top left: 6*0 = 0.1. Top right: 
60 = 0.03. Bottom left: 60 = 0.01. Bottom right: 60 = 0.003. 
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Fig. 31. — e!^^-{D) ( "Ghirlanda" ) plots for universal jet models based on top-hat profiles. 
Each panel shows the result for 7-^ = 0.1, 0.03, 0.01, and 0.003. Top left: 60 = 0.1. Top 
right: 60 = 0.03. Bottom left: 60 = 0.01. Bottom right: ^0 = 0.003. 
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Fig. 32. — dN / dlog{Eiso) distributions for universal jet models based on top-hat profiles. 
Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: ^0 = 0.1. Top right: 
6*0 — 0.03. Bottom left: 9o — 0.01. Bottom right: 9o — 0.003. The vertical scale is arbitrary. 
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Fig. 33. — dN/d\og{Ej'') distributions for universal jet models based on top-hat profiles. 
Each panel shows the result for 7 = 10, 33.3, 100, and 333. Top left: 6q = 0.1. Top right: 
9q — 0.03. Bottom left: 9q — 0.01. Bottom right: ^0 — 0.003. The vertical scale is arbitrary. 



